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With the developments of large galaxy surveys or cosmic shear surveys it is now possible to map 
the dark matter distribution at truly cosmological scales. Detailed examinations of the statistical 
properties of the dark matter distribution reveal the detail of the large-scale structure growth of the 
- - , Universe. In particular it is shown here that the behavior of the density field bi-spectrum is sensitive 

\l . to departure from normal gravity in a way which depends only weakly on the background evolution. 

' The cosmic shear bispectrum appears to be particularly sensitive to changes in the Poisson equation: 

we show that the current cosmic shear data can already be used to infer constraints on the scale of 
a possible higher-dimensional gravity, above 2 ft~^Mpc. 

^ • PACS numbers: 98.80.Cq, 98.80.Es, 04.80.Cc, 04.50.-fh 
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Almost all cosmological observations point toward a coherent model in which there is a non- vanishing dark energy Q . 
■ The existence of such a component is very puzzling from a high energy physics point of view and signals physics beyond 
fSl the standard model in way we do not yet apprehend. The resolution of such an issue may lie either in the introduction 
of a new form of cosmic fluid - a dark energy - with a negative pressure that would correspond to two third of the 
actual energy density of the universe or maybe in a modification of the Einstein gravity equations at cosmological 
^-p • scales. The two scenarios have obviously to be distinguished 0]! The first route is usually taken in models evoking 
quintessence component. However it cannot excluded that a non-zero cosmological constant observation actually 
signals a modification of the Einstein law at cosmological scales. This idea is actually fuelled by recent developments 
in string theory phenomenologies. Indeed the introduction of branes in the context of superstring theories0,|3|, led to 
concepts of higher dimensional spacetimes. Usually it is assumed that the interaction gauge fields are localized on a 
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O ■ 3-brane (i.e. a 3-1-1 dimension spacetime) whereas gravity propagates in all dimensions. In any of such string inspired 
^ , models, one expects both the existence of Kaluza-Klein gravitons implying a non standard gravity on very small scales 
c/2 ■ and light bosons, which can manifest as a new fundamental small scale force. It cannot be however excluded that 
^ ' gravity might be modified at large cosmological scale as well 0, IT^ . This possibility has been explicitly linked 
to the evidences of a cosmological constant J5j. It asks for precise test of gravity law at cosmological scales. This 
is the aim of this letter in which we pursue investigations on how cosmological observations, namely cosmic shear 
observations, could be used to test gravity on large distances. 

Deviations form Newtonian or Einstein gravity have been looked for in the past mainly at rather small scales 
(from a cosmological point of view) whether it is in laboratory experiments 0) 12] or in stellar systems. That put 
severe constraints on the post- Newtonian parameters 0, ^3 but tells little on the validity of the Einstein equations at 
cosmological scales (say at scales above 1 h^^Mpc). At best one can invoke studies of galaxy clusters. Comparisons 
between X-ray emissivity and gravitational lensing, which appears to be an indirect test of the Newton law through 
the equation of hydrostatic equilibrium, show no dramatic discrepancy below 2 Mpc if normal gravity is assumed. 
On larger scales, there is no possible test on the Poisson equation (although other aspects of GR can be tested, 
see '16]) but by the mechanism of structure formation through the development of gravitational instabilities. This is 
the object of this letter. 

Following we simply assume that the gravitational force is changed at large-scale. Other possible modifications 
of the large-scale gravity have been considered in W§\ . We assume in particular that the background metric is that 
of a Friedmann-Lemaitre metric with an expansion factor o(z) having a redshift dependence assumed to be that of a 
Univers with a cosmological constant. Note that it does not mean the Friedmann equations is assumed to be correct, 
e.g. that there is a dark energy component, but that we take the evidences for a non-zero cosmological constant at 
face value. The latter being mainly of cosmographical nature, e.g. they come from the redshift dependence of the 
expansion parameter and quantities that are directly related to it such as angular or luminosity distances, they do 
not provide us with an actual weighing of the matter/energy content of the Universe. 

As a further consequence, in what follows, as long as we are dealing with subhorizon scales, we can take the metric 
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to be of the form, 

ds^ = -(1 ~ 2$)di2 + a2(l + 2$) (dx' + g'(x)df^') , (1) 

where t is the cosmic time, a{t) the scale factor, x the comoving radial coordinate, dfl^ the miit solid angle and 
q{x) = (sinx, X, sinhx) according to the curvature of the spatial sections. In a Newtonian theory of gravity, $ is the 
Newtonian potential <i> n determined by the Poisson equation 

A$Ar = inGpa^S (2) 

where G is the Newton constant and A the three dimensional Laplacian in comoving coordinates, p the background 
energy density and S = Sp/p is the density contrast. If the Newton law is violated above a given scale then we 
have to change Eq. and the force between two masses distant of r derives from $ = $Ar/(r/rs) where f{x) 1 
when X ^ 1. This encompasses for instance the potential considered in 0,123 for which f{x) ~ 1/(1 + x) (in that 
case f (xl/x and 5D gravity is recovered at large distance). Using Q it leads, with r = ax, to 

which, making use of A[f{x)/x] — — 47r(5(^^(x) + fs{x/xs) with fs{x/xs) = {d^f )/x, gives 

A$ = A$Ar - Gpa^ J d^ii'S{x' + x)fs{x'/xs). (4) 

The motion equations of the field are better handled in Fourier space. In Fourier space, Eq. ^ reads 

- fc2$(fc) = 4TrGpa%k)fc{kr,) (5) 

where /^(fcr^) = 1 — 27r^ fs{krs), fs{krs) being the Fourier transform of fsir/rg). It is unity for large values of fcr^, 
it behaves like k for small. 

To explore the physics of large-scale structure growth the now standard method is to expand the local density 
field with respect to the initial density field. Such an expansion can be made either in real space or in Fourier space, 

5(k,t) =5(i)(k,t)+(5(2)(k,t) + ... (6) 

where (5'^^)(k) is the linear density field, ^'^^^(k) is quadratic in the linear density field, etc. The motion equation for 
the density field can be solved order by order and the resulting expansion can be used to infer the evolution of the 
stochastic properties of the local scalar field, whether it is the local density field or the gravitational potential. More 
quantitatively for any stochastic field X (statistically homogeneous and isotropic) we define its power spectrum Px 
by 

(X(k)X(k')) =7'x(fc)<5(3)(k + k') (7) 

where J^^-* is the Dirac distribution, X the Fourier transform of X and the brackets refer to ensemble averages 
over the stochastic process that gave birth to X. Note that the poisson equation then implies, 

VA^^{k)^{47:Gpa^frs{k). (8) 

It is usually assumed that the initial metric fluctuation are Gaussian distributed. In the early stage of the dynamics 
it is therefore enough to consider the power spectrum for the scalar fields to fully characterize them. The subsequent 
evolution of the field however is bound to induce non-Gaussian effects. The second order 5^"^^ (k, t) in particular 
contains mode coupling terms that induce non-zero third moments (whether it is in real or in Fourier space). One 
way to quantify those effects is to introduce the field bispectrum, i3x(ki, k2, ks) defined as, 

(X(ki)X(k2)X(k3)) =Sx(ki,k2,k3)<5(3)(ki+k2+k3). (9) 

In cosmology it is common to define the reduced bispectrum as. 
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the reason being that turns out to be roughly scale independent. 
With modified gravity laws we can already note that, 

V^^{k) = {^TzGpa'fvmfcikr.f (11) 

and 

for equilateral wave vector configurations. 

As a consequence, a simple way to test the validity of the Newton law is to compare power spectra of the density 
field and that of A<I> to see if they follow the same scale dependence. That was the point made in This is 

possible however only if one can measure 5 and <& independently. That supposes in particular that galaxy counts can 
give a reliable account of the mass distribution at large-scale, which as some level is always questionable due to the 
unavoidable presence of biasing effects 1311 Note that constraints on modified gravity law can also be obtained from 
the shape of the power spectrum alone |2q | but such constraints rely on assumptions on the shape of the spectrum of 
the primordial metric fluctuations. 

The aim of this letter is to show that it is possible to test departure from standard gravity from cosmic shear survey 
alone taking advantage of the behavior of the reduced bispectrum. 

As long as we have not entered too much into the nonlinear regime, the amplitude and shape of the power spectrum 
is determined by the growth rate of the linear term, (J*^^-* (k, i). In such model, and contrary to cosmological models 
with standard gravity, the growth rate is scale dependence. More precisely we have, 

5^-^\\i,t) ^ Dkm^.O^) (13) 

where Dk (t) is the linear growing rate of mode k that satisfies the equation j23| , 

bk{t) + ^Hbkit) = ^H'' n{t) /, {kr,/a{t)) Dk{t). (14) 

The solution of this equation have already been explored in • It leads to a slowing of the growth rate for the scales 
beyond Tc. 

For the 2nd order terms the structure of the conservation and motion equations leads to the following functional 
form (see also the approach adopted in |27||). 



5(2)(k) = |d3ki|d^k2^2i£^^L±M[F„(ki,k2)4n.(ki)<5n„.(k2)- 



Fi(ki,k2)<5lin.(ki)<5lin.(k2)^^ +F2(ki,k2)4n.(kl)<5li„.(k2) ^ J 



(15) 



where the coefficients Fq, Fi and F2 represent the amplitude of respectively the monopole, dipole and quadrupole 
terms that appear of the expression of the 2nd order local density field. For cosmological models with a standard 
gravity, Fq, Fi and F2 are scale independent. Moreover they appear to behave like D^{t) with a fixed coefficient that 
depends extremely weakly on the cosmological parameters. 

Lest us now explore the case of a modified gravity dynamics. The continuity equation imposes that (^'■^■') = 
which implies that Fq — Fi + 2/3F2 = 0. The behavior of the 2nd order local density field will then be obtained from 
the evolution equation for these 2 functions. 



Fo-2HFo-^H^mfc[\^i+'i^2\^^]F„ = -b^^b^^ + -H 



-a{t) 



■a{t) 



Dk.Dk, (16) 
(17) 



and 4/3£'fcjI?fc2 in the r.h.s. of the previous equation is changed into 2Dk-^^Dk^ for the evolution equation of Fi. 

These two equations have a very similar structure and can be solved numerically for whichever background we 
consider. To get the behavior of Fq and Fi in the following we assume we live in a universe in which, 



a{t) — ao sinh 



(18) 
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FIG. 1: Time dependence of Fo(ki,\i.2) / D^-^ / Dk^ for fci = k2 and for different value of the relative angle between ki and k2. 
This angle varies from to tt from top to bottom lines. 
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FIG. 2: Shape of the bispectrum for the density field (long dashed lines) and the potential Lagrangian (solid lines). The 
computation has been made for fci = ra/10 and k2 — 2rs/10 and assuming P{k) ~ k^^ for the density field. 

which corresponds to the expansion factor behavior for a universe for a given value of fio and no curvature (in the 
following we do the calculations with = 0.3). Note that we found that the results described in the following are 
very weakly dependent on the background behavior. 

Compared to the standard gravity case, _Fo, ^^i and F2 now depend on the wave vectors through both their length 
and their relative angle. It introduces further dependence of the bispectrum on the triangle configuration. As an 
illustration we present on Fig. ^the behavior of fo(ki, k2)/-Dfci/-Dfe2 as a function of time and for different angle 
between the two wavevectors (of the same length). The angular modulation is exhibited at late time (at early time 
this ratio is 17/21 irrespectively of the angle). It is found to be actually modest given the precision with which such 
a quantity can be measured. Similar results can be found for Fi . 

The resulting behavior of the bispectrum are show on Fig. |5] It shows the time evolution of the reduced bispectrum 
Q as defined in Eq. (|10|l for either the density field or the potential Laplacian compared to the standard gravity case. 
The density bispectrum can be measured for instance in galaxy catalogues. It has been done in particular in the 
PSCz galaxy survey [2^. It does not show however strong dependence on the modified gravity efi^ects. This is not 
the case for the Laplacian field whose reduced skewness is amplified by modified gravity effects. That was indeed to 
be expected from the scaling relation shown in Eq. (|12|l which suggests that such a quantity is expected to grow at 
large scale. 

It happens that quantities directly related to the reduced third moment of the potential Laplacian have been 
measured in cosmic shear catalogues l2^. The three-point correlation of the shear field (or rather some geometrical 
average of it) has been detected with a significant confidence level (through a rather elaborate procedure) at different 
angular scales [s^l- If one applies the scaling law one expects for modified gravity to the predictions of standard 
gravity (in a A dominated universe) it is possible to put constraints on Vg. From Fig. (Sjit appears that current data 
already disfavor values of below 2h~^ Mpc although it is probably impossible to put a definitive statement from 
this data set. The constraints on the scale at which modified gravity intervenes are however strongly dependent on 
the modelling of the observations, and on the detailed shape of the 4D to 5D transition. 

Cosmic shear surveys appear nonetheless very attractive to test the validity of the Einstein gravity law at cosmo- 
logical scales. The fact such surveys provide us with genuine potential distribution at large scale makes them precious 
probes of the gravitational instability mechanisms in general. Undoubtedly, much stringent constraints should be 
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FIG. 3: The measured reduced shear skewness (square with error bars) compared to cosmological models with standard gravity 
(dotted line for A-CDM, upper dashed line for open CDM and lower dashed line for standard CDM) or with modified gravity 
for A CDM, Ts = 10/i~^ Mpc (lower solid line) and Ts — 2h~^ Mpc (upper solid line). 



obtainable from the coming cosmic shear surveys such as the CFHTLS or on the long term satellite missions like 
JDEM/SNAP. 

Note that should such a modified gravity be the reason of the observed evidences of an accelerated universe would 
be expected to be the order of the Hubble size. Obtaining constraints on of this order form large-scale surveys may 
not be as hopeless at it seems. For the potential shape derived in 14] parameterized by a scale Tc whose estimated 
value is Vc — 1.2 c/Hq 0|, the amplification of the three-point correlation is not as small for the available scales as 
the simpHfied model used here suggests. At 10 h~^Mpc scale the effect would be at percent level; at 100 h~^Mpc 
scale, it would reach 7% (a valid theory of gravity has however to be done in those models). If this approach reveals 
the only way to distinguish between quintessence type and modified gravity models to account for the dark energy 
evidences, it calls for the realization of very large - possibly whole sky - cosmic shear surveys. 
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